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Abstract 

The paper addresses the problem of a a semi-infinite plane crack along the interface between two 
3D isotropic half-spaces. Two methods of solution have been considered in the past: Lazarus and 
Leblond (1998) applied the "special" method by Bueckner (1987) and found the expression of the 
variation of the stress intensity factors for a wavy crack without solving the complete elasticity 
problem; their solution is expressed in terms of the physical variables, and it involves five constants 
whose analytical representation was unknown; on the other hand the "general" solution to the 
problem has been recently addressed by Bercial-Velez et al. (2005), using a Wiener-Hopf analysis 
and singular asymptotics near the crack front. 

The main goal of the present paper is to complete the solution to the problem by providing 
the connection between the two methods. This is done by constructing an integral representation 
for the Lazarus-Leblond's weight functions and by deriving the closed form representations of the 
Lazarus-Leblond's constants. 
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1 Introduction. 



The fundamentals of the theory of three-dimensional interfacial cracks and the corresponding integral 
equation formulations were introduced by Willis (1971a, 1971b). A singularly perturbed problem of a 
general static loading of a semi-infinite plane crack along the interface between two isotropic half-spaces 
of different linear elastic materials is considered here. The two characteristic features of this problem 
are the coupling of all three opening modes and the oscillatory behaviour of the solution near the crack 
edge. For this problem two issues are of major interest: the representation of the stress intensity factors, 
and the expression of the variation of stress intensity factors arising from an infinitesimal coplanar 
perturbation of the crack front. Two methods of solution have been considered in the past: Lazarus 
and Leblond (1998a, 1998b) applied the special method by Bueckner (1987) and found the expression 
of the variation of the stress intensity factors without solving the complete elasticity problem; their 
perturbation formulae are simple and elegant, expressed in terms of the physical variables, but these 
formulae involve five constants whose analytical representation was left unknown; on the other hand the 
general solution to the problem has been independently addressed by Antipov (1999) and Bercial-Velez 
et al. (2005). 

In the present paper, we revisit the work by Bercial-Velez et al. (2005) and show that it is possible to 
construct the weight functions in the form directly related to the weight functions of Lazarus and Leblond 
(1998a). This also yields the required analytical expressions for the unknown constants in asymptotic 
formulae of Lazarus and Leblond (1998a). Then, using an asymptotic method and the integral reciprocal 
identity, formerly introduced by Willis and Movchan (1995), we obtain canonical integral representations 
for the stress intensity factors. 

The structure of the paper is as follows. Section 2 includes the governing equations and presents 
the main result. The fundamental identity and the weight functions are described in Section 3. Closed 
form representations for the Fourier transforms of the weight functions are given in Section 4. Local 
asymptotics and comparison with the Lazarus-Leblond's weight functions are given in Section 5. Finally, 
Section 6 includes a brief outline of applications to the wavy crack problem. 

2 Problem definition and main result. 

We consider an infinite elastic body, consisting of two different isotropic materials that occupy the 
upper half-space (material 1) Rg = {(xi, x 2 , x 3 ) £ R 3 : x 2 > 0} and the lower half-space (material 2) 
Rg" = {(xi, x 2 , x 3 ) S R3 : X2 < 0} (see Fig. [T]). The Poisson's ratio and the shear modulus of materials 
1 and 2 are denoted by 1/+, fi + and fi-, respectively. 

The crack lies on the half-plane R^ = {{xx, £2, £3) £ R3 : x\ < 0, X2 — 0}, whereas the ideal contact 
conditions are valid on the half-plane R^ = {(x\,X2.,x 3 ) € R3 : x\ > 0, X2 — 0}. This implies that the 
traction and displacement components are continuous across the interface ahead of the crack front 

[0-12] = [0-22] = [0-32] = 0, [ui] = [u 2 ] = [u 3 ] = 0, as xi > 0, (1) 

where the square brackets denote the jump while crossing the interface, 

[f]{xi,x 3 ) = f(x 1} +0, 13) - f(x u -0, x 3 ). (2) 

The Stress Intensity Factors Kj, j — I, II, III, are defined as in Hutchinson et al. (1987) (also see Lazarus 
and Leblond (1998)), so that the asymptotics of tractions ahead of the crack edge and the asymptotics 
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of the displacement jump across the crack surfaces take the form 



K(x 3 ) 

CT22 {xi, 0,x 3 ) +icr 12 (x 1 ,0,x 3 ) ~ x\ e , xi 

V Zttxi 

, n \ ^III (£3) n+ 



o+, 



(3) 



and 



V27TX1 



[U2 + IU 1 }{XI,X 3 ) /o ; — K{x 3 ) 



(1/2 + ie) cosh(7re) 



-xi 
2ir 



x\ -> , 



[u 3 ] (2=1,^3) ~ 2 ( — + — ) ^iii(x 3 ) 



(4) 



-x\ 
2ir ' 



a;i -> 0" 



respectively, where if (x 3 ) = ifi(x 3 ) + iifnfe) is the complex stress intensity factor and 

1 „ M+ + (3 - 4^+)/i- 



6 2tt l0g /x_ + (3 - 4i/_)/i+ 



(5) 



is the bimaterial constant. 

The crack faces are loaded by surface tractions with components — Pi(x\,x 3 ) on the upper face 
(x2 = + ) and Pi(x\,x 3 ) on the lower face {xi = CP). 

The body forces are assumed to be zero, and the displacement components satisfy the Lame system 



u k ,ik + (1 - 2u±)u. 



i,33 



0, in R± 



(6) 



with boundary conditions 0^2(2:1, X2, x 3 ) — Pi(x\,x 3 ) as x 2 — > O*, on the crack surfaces (21 < 0), and 
ideal contact conditions [uj](a;i, £3) = [042] (£1, £3) = on the interface (x\ > 0). 
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Lazarus and Leblond (1998) derived expressions for the variation of SIFs generated by a small copla- 
nar perturbation of the crack front, x\ = 5<j){x^). In terms of the Fourier transforms and assuming that 
K and Km are initially uniform, Lazarus-Leblond's formulae are written as follows, 

A*(A) = - 1+*L ( 7+ 5^«r(l - 2ie)^^ + ^K\X\ 

y ' da v 8cosh(7re) \ ,+ e V ' 1 + 2ie ' 11 

47m cosh(7rc/2) . , JA| 1+ * ) ~ 

r(l - ze Xmsign A — — — f (50, (7 

1 — 1/ e 1 + ze 



... d-Kin . , ( Try T , ,,, . 1 — ^cosh(7re/2) . 

Aif m (A) = — T— <50 - ^ —p -Kin A - ^ ^ sign A Im 

da [_ 4 2 e 



1 + ie 



where the tilde denotes Fourier transform with respect to X3, the star denotes complex conjugation, 
dK/da and dKm/da represent the derivative of K and Km, respectively, for a uniform advance of the 
crack front and Y is the gamma function. 

Here 7+, 7 _, 7m, j z are complex constants, and 7 is a real constant, and the analytical representa- 
tions for these constants were absent in the literature. We note that the constants 7+, 7_, 7m, j z , 7 
are taken from the asymptotics of the weight functions used by Lazarus and Leblond (1998). 

The main aim of the present paper is to fill-in this gap and to obtain the representation for the 
Lazarus-Leblond's constants. This was possible due to the analytical results of Antipov (1999) and 
Bercial-Velez et al. (2005). The required formulae have the form 



2 3(6 + e) - x/o 2 - d 2 

7 ~ 



7r ^Jb 2 - d 2 + b + 



e 



7111 



8^e(l + ie)Vb 2 -d 2 Vb 2 -d 2 



2*(1 + 2te)r(l/2 + ie)r(l - ie)(V6 2 - d 2 + b + e)(Vb + d + Vb~d) ' 

lz = - 7in{1+ 2ie)-^=, (8) 
86 

7- 



7+ = 



tt(1 + 2ie)( v / 6 2 - d 2 + 6 + e) ' 
4e6r(l/2 - ie)(Vb 2 - d 2 - 6 - e) 



where 



4* e aT(l/2 + ie)r(l - 2ie)(Vb 2 - d 2 + b + e) 
1 — Va_ 1 — v_ , 1 — 2z/_i_ 1 — 2z/_ ^ 1 i/_ 



/Lt+ /Lt_ 2/X+ 2/i_ /Lt+ jU_ 

Note that, using the parameters in ([5]), the bimaterial constant e given by © can be written as 



(9) 



1 1 b + d nn\ 
2^ l ° g —d- (10) 



1 These expressions were kindly supplied to us by J.B. Leblond. 
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3 The fundamental identity. 



We will adapt the method introduced for the homogeneous elastic space by Willis and Movchan (1995). 
This method involves the use of the reciprocal theorem (Betti formula) in order to relate the physical 
solution to the weight function, which is a special singular solution to the homogeneous problem. 
In the absence of body forces, the Betti formula takes the form 



{al>n j u?>-al>n j u?>}ds = 0, (11) 
on 



where d£l is any surface enclosing a region f2 within which both u\ and u\ satisfy the equations of 

(1) A (2) 



equilibrium ©, with corresponding stress states and <J^\ and rii denotes the outward normal to 



on. 

Applying the Betti formula to a hemispherical domain in the upper half-space R^~, whose plane 
boundary is x 2 = + and whose radius R will be allowed to tend to infinity, we obtain, in the limit 
R — ► oo, 

J {ag 5 + > x 3 )u[ 2) {x u 0+, x 3 ) - 4 2) (zi, + , x 3 )uW {xuO+^s^dxtdxa = 0, (12) 
(2:2=0+) 



provided that the fields u\ and u\ 2 ' decay suitably fast at infinity. We can also assume that 
represents a non-trivial solution of the homogeneous problem, whereas 16 stands for the physical field 
associated with the crack loaded at its surface. 

Similar to Bercial-Velez et al. (2005), we now define a new vector function {CA}f =1 in the following 
way 

U 1 (-x 1 ,x 2 ,-x 3 ) = -u\ '(x 1 ,x 2 ,x 3 ), 

Ua(-xi,X2,-x 3 ) = u 2 2 \xi,x 2 ,x 3 ) 1 (13) 

Ua(-xi,X2,-x 3 ) = -u 3 > (x 1 ,x 2 ,x 3 ), 

which corresponds to introducing a change of coordinates in the solution u' 2 ' , namely a rotation^] about 
the X2-axis through an angle w. It is straightforward to verify that the function Ui satisfies the equations 
of equilibrium ([S]), but in a different domain (see Fig. In the sequel, the vector function {[7j}f_ 1 will 
play the role of the weight function, whereas the vector function {it[ 1 ''}f =1 will be identified with the 
physical solution (and we will drop the superscript (1), no longer needed in the following notations). 
The notation Y^hk(Xi, X 2 , X 3 ) will be used for components of stress corresponding to the displacement 
field Ui(Xi, X 2 , X 3 ). Equivalently (see (fT3)) ), the components of displacement and stress are related by 
the formulae 

u\ '(x 1 ,x 2 ,x 3 ) = R lh U h (-x 1 ,x 2 ,-x 3 ), (14) 

and 

(21 

al j '(x 1 ,x 2 ,x 3 ) = R ih 'Shk(-xi,X2,-x 3 )Rkj, (15) 

where 

-10 

R= ( 1 1 . (16) 
0-1 



2 It is noted that this transformation of coordinates differs from the transformation of reflection used by Willis and 
Movchan (1995) for the case of a crack in a homogeneous elastic space. 
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Replacing u\ (xi, X2, X3) with u\ (xi — x'i,X2,X3 — x 3 ), which corresponds to a shift within the plane 
(xi, £3), we obtain 

Wi2(xi,0 + ,x 3 )Ri h Uh{x' 1 - xi,0 + ,x' 3 - x 3 ) 

-Ri h T, h 2{x[ - a;i,0 + ,X3 - x 3 )u i (x 1 ,0 + ,x 3 )}dx 1 dx 3 = 0. (17) 

A similar equation can be derived by applying the Betti formula to a hemispherical domain in the 
lower half-space Rg", 



{(T i2 {x 1 ,Q' ,x 3 )R ih U h {x' 1 - xi,0~,x' 3 - x 3 ) 

02=0-) 

-~R ih ^h2(x' 1 - xi,0~,x' 3 - x 3 )ui(xi,0~ ,x 3 )}dxidx 3 = 0, (18) 
and hence by subtraction of (fTS)) from (JT7J) , we obtain 



{a i2 (xi,0, x 3 )R lh [U h ](x' 1 - xi,x' 3 - x 3 ) 

-Rih^^ix^ - xi,0,x' 3 - xz){v.i](x],x 3 )}dxidxz = 0, (19) 
which is a convolution integral and the equivalent representation is 

{<Ji2 * Rih[U h ] - fii/,S/,2 * Wj^i.ia) = 0. (20) 

Upon writing 

a l2 {xi, Q,x 3 ) = pi(x!,x 3 ) + a { ^\xi,x 3 ), (21) 
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where Pi(xi,xs) = tr^^i, 0, X3)H(— x\) is the loading, and cr> 2 '{x\, X3) — ai2{xx,Q,x 3 )H{xx) is the 
stress field ahead of the crack edge, with H(x\) being the Heaviside function, we get the Betti identity 
in the following form 

K ( 2 +) * Rihph] - Ri h V h2 * MIW, 4) = -{Pi * RihWhW^x's). (22) 

Equation (|22[) will be used for evaluation of the stress intensity factors. A field with components 
Ui(x\, X2, X3) is the weight function defined as follows 

(a) it satisfies the equation of equilibrium ([6|); 

(b) [Ui] = when x x < 0; 

(c) the associated is continuous and £^2 = when xi = and x\ > (homogeneous boundary 
conditions); 

(d) [Ui] ~ k i (xi)x 1 i ^ 2 5(xa) as x\ —* + , where ki(x%) is a bounded function, and 5{x 3 ) denotes the 
Dirac delta function|f| 

(e) Ui is a linear combination of homogeneous functions of degree —3/2, —3/2 + it and —3/2 — ie. 

4 The weight functions. 

Let us introduce the Fourier transforms of the displacement jump and of the traction components for 
the weight functions 

/OC /"OC 
/ [U i ](x 1 ,x 3 )e i ^+ iXx3 dx 1 dx 3) 



S i2 (/?,7) = 




The weight functions are defined in such a way that [Ui] + is analytic in the upper half-plane 

C+ = {-00 < Re(/3) < 00, Im(/3) > 0}, (24) 

whereas T, i2 is analytic in the lower half-plane 

= {-00 < Re(/3) < 00, Im(/3) < 0}. (25) 

4.1 The Wiener-Hopf problem. 

The relationship between the Fourier transforms of displacement jump and traction components has 
been derived by Willis (1971a, 1971b) and takes the form 

[U)+(P,\) = ±G(P,X)X-(P,\), (26) 
P 

3 The bounded functions ki(x\) have the form ki{x\ ) = k^xf + k^x^ te + fc' 3 ' , where k^\ fc' 2 ' , fc' 3 ' are constants. 
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where A 6 R and [i7] + (/3, A), X (/3, A) are the limit values of the functions as Im(/3) 

/ -idf3p bp 2 + e\ 2 -e(3X \ 

bp 2 id(3p 



p=^/3 2 + X 2 , G 



idXp 

V -idXp -e/3X bp 2 + e/3 2 j 



l-u+ 1 - v- , l-2v+ 1 - 2z/_ 
b = 1 , d = 



M- 



2p, 



2^_ 



Without loss of generality we assume d > 0. 

Note that the order of the components of the traction vector £ in equation ([26 
the standard one, namely 



^22' ^12 J ^32 



Following Antipov (1999), the matrix G can be written as G = J1G0J2, where 



1 



1 



G (1 



Ji = J^ T =(100, J 2 = J2 1 = 1 

010/ v 1 

/ bp 2 idXp id(3p \ 

-idAp fop 2 + e/? 2 -e/3A 

V -idf3p -ej3X bp 2 + eX 2 J 
Here we use the following normalization with respect to mechanical parameters 

4, + (P 1 X) = J^ 1 [U]+(p,X), c/>-(P,X) = bJ 2 W(0,X), 
so that the problem can now be written as 



where 



0+(/3,A) = -G 1 ( / 9 ! A)0-(AA), 



/ p 2 id*Xp id*{3p \ 

—id*Xp p 2 + e*/3 2 — e*/3A 



\ -id^fip -e*/3A p 2 + e*A 2 / 
and the dimensionless parameters d* and e* are defined as follow: 

j d e 
a* = - < 1, e* = -. 
b b 

Using a new variable £ = /3/|A|, we can write equation (I3ip in the form 



^+(0 = -G.(e,Bign(A))^(0, 



where 



p* = + 



pi 



id* sign(A)p* id*(,p* \ 



-id* sign(A)p* p 2 + e*£ 2 -e*£ sign(A) 



V —id*£p* 



-e*£sign(A) 



pl + e* 



(35) 



4>t(0 = + (C|A|, A), - — 0-(e|A|,A). 



Note that is analytic in Ct , whereas is analytic in C f . The solution of the Wiener-Hopf equation 



(j34|) is outlined in the Appendix. 

4.2 Fourier transforms of the weight functions. 

The Fourier transforms of three linearly independent weight vector functions and the corresponding 
traction components are: 

(1) the first weight function 



ti = [ 




+ 








£ = [ 




■ 








# = [ 




■ _ 









-1/2 



+ 2F x (0f (€ + i)- 1/2 , 



|A|' 



sign(A) 



(l + e.)(€-0 



+ 2td.ei\(0f(C-i) 1/2 . 



mSi = -2*d, sign(A)P 2 (0(£ - *) 1/2 , 



As- 1 - 



(l + e.)(£-i) 



2id.iM0}(£-i) 1/a ) 



(2) the second weight function 



*3 




i a = 


\ul 2 ] + = 


t 2 2 = 


[ul 3 } + = 


= 


iAr* 12 


fc 2 = 




^2 = 


As- 2 

1 i ^*32 



1 + e* 



+ ^ 2 (0 (£ + i)- 3/2 , 



sign(A) 



1 + 



2dgd*(i — a) 



1 



e-^ 2 (oj(? + *)- 3/2 , 
(e-*) 1/2 



2dgd*(i - a) 



-rf*p 2 (0^i(0 



(£-*) 1/2 



(e-*) 1/2 



2 = ^ Sign(A) 1 2dgd.fr -a) + d * P ^\ £ + 



(36) 



(37) 



(38) 



(39) 
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(3) the third weight function 



'*3 



d*(£, - a) 



(40) 



* 3 |A| * 12 \ £-<z <z) i.(t-o). 

t - 3 _ » E -, _ |m» + ftKL _ J (£ _ if n (41) 

|A| I £ ^ a — a) a*(2 — a) 



,- 3 & v" 3 • • rx^ / ^ 2 (0 _l jMO 

0*2 = TTT S *32 = sl g n ( A \ — 7. + -TTc V ~ 77 

|A| L ? — a «*(£ — a) a*(z — c 

Here, the functions Fj,Pj, j = 1,2, are defined by 

; M- A+(QS+(0 .M_ + 

Pl (0 = M- A * (0 7f*" (0 , P 2 (0 = M_A-(OcosB-(£). 

MO 



(42) 



The functions , , Aj and the constants a, M_, <io are defined in the Appendix, equations (|A.5|) 9. 
(TO) , and (TA~7l) . 

5 Comparison with the Lazarus-Leblond weight functions. 
5.1 Local asymptotics. 



Let us consider the Fourier transform of the fundamental identity 

<ft 2 Rih\ui}+ - flihE^fo]- = -ft^ h [C/i] + . (43) 

Using the asymptotics for the physical fields and weight functions given in the Appendix, we can 
observe that the structure of the asymptotics for £ is as follows 

a+ = A a p-^ 2 + A+fr 1 /^ + A-p- 1 ' 2 -™ + o(|/?r 3 / 2 ), 

- co, (44) 

[U]+ = A [v] 0- V» + A+ v] 0-y^ + A- m 0-y^ + 0(|/3|- 3 / 2 ), 

where A a ,A^ and A[u] , -Af™ are complex quantities depending on A. Note that in equations (|44p and 
in equations (|47p in the text below we suppress the subscript and superscript indices (compare with the 
identity (|43| ) to simplify notations. It follows that 

^ 2 RihPi\ + ~ A jk (\)K k (\)j^, /3-oo, 0£C+, (45) 
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where {K k (X),k — 1,2,3} = {K(X), K* (X), Kiu(X)}, with the tilde denoting Fourier transform with 
respect to x% only, Ajk(X) = y / |A[{a J fc(A)}/(4v / 2) and 

... . ^d.VT^Io D ... . ... 2d*VT=ia~\X\ ie 

an(A) = -sign(A) — -, 012(A) = -sign(A)- 



ci |A|"' w & w c 2 £> 



.,. d*\/l — ia D d*-*/! — * a |A| l 

021 (A) = 777777 022(A) - 



ci |A| i6 ' " w c 2 D 1 

2(1 -ia)dl' 2 D 2(1 - iaJdyVl - ad* IAI* 

a 3l(A) = - ; ., /, , TTTTTi a 32 (A) = 



(46) 



(l + i)ciVl-ad» |A| ie ' " w (1 - i)c 2 D ' 

a 13 (A) = 2^(1 + 0, a 23 (A) = - sign(A) ^l±lKl + gll a 33 (A) = 0. 

We also note that, for /3 S , 

£~ = As/3 1 / 2 + A+/3 1 /2+ i e + 4-/31/2-* + 0{\p\- l l 2 ), 

/3^oo, (47) 

[«]- = A N r 3/2 + 4]/?- 3/2+ie + ^ M r 3/2 - 4£ + oa/?r 5 / 2 ), 

and hence 

fl^^r-A-tW^iA)-^, /3-»oo, /JeC r (48) 
Then similar to Willis and Movchan (1995) we can write 



at 2 Rih[Ui} + - RihTi^Ui}- ^A jk {X)K k {X)\^-^-j— -1, /3^oc, (49) 

where the term in brackets can be identified as the regularization of the delta function, namely —2i:iS(P). 
Equation (|4"3")) implies 

lim ^[a+RihpiV - RihXi'iui}-} = -iA jk (X)K k (X) = - lim ^[p^p 3 ^], (50) 

x\— >0 1 ij- >0 1 

and hence the Fourier transforms of the stress intensity factors become 

r° 

K k {X) = -i lim A^(X) / p i {x 1 ,X)R ih [Ul}+(x , 1 - x 1 ,X)dx 1 , (51) 
^i^o j_ oa 

where Ajy(A) = B/y(A) = \/2{6fcy(A)}/ [Vl - ia(-ad* + 1 + e*)] and 

611(A) = -sign(A) ■} — — — -==, 012(A) = 2cx[-l + {-1 + a)d*\a- 



1 1 



621(A) = -sign(A) — — - — —2 — — 7==, o 22 (A) = 2c 2 [l + (-z + a)oV~ 



{-i + a)dl |A|«vW |A|"vlAf' 



631(A) = -ii-^ad*^^, 6 32 (A) - -V2sign(A)(l +z)(l + za)ad 2 ^^, 
v2 V|A| Vl A l 

6i 3 (A) = 2ici[(l + m)d 2 + e*]^=, 623(A) = 2ic2[(l + m)^ -2-e,] , (52) 

£V|A| |ArV|A| 



633(A) = -\/2sign(A)(l + i)(-< + "H* ^^ ' 

V |A| 
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Note that the components of the matrix B are given by B k j (2:3) = -\/2{ 6^(2:3)} / [Vl — ia(— ad* + 1 + e*)l 
where 

, , N i . , ,ci[-l + (-i + a)d,](l + e*)cos[7r(3 + 2ie)/4]r(l/2 + ie) 

611 2:3 = sign 23 — , 

7T (-« + a)^ Dlzal^Vksl 

, , , 2 . 1 . . ... cos[7r(l + 2ie)/4]r(l/2 + ie) 

, / s * • / .c 2 [l + (-« + a)cZ*](l + e*) 1 ■ w^itvi /o • \ D \ x 3\ le 

621(2:3) = sign(x 3 ) r — - — cos tt(3 - 2ze)/4 r(l/2 - ie) 

7r (—i + ajaj 



T 3 



622(2:3) = -02 [1 + f-i + a)d,}a cos[tt(1 - 2ie)/4]r(l/2 - 



631(2:3) 



1 — i yl — ia 
ma* — , , 



20? 



Fa 



2z 



613(2:3) = — ci[(l + ia)rf» + e s 

7T 



632(2:3) = — p= sign(x 3 )(l + i)(l + ia)aa\,, 

v 7 *" VFsl 

cos[tt(1 + 2ie)/4]r(l/2 + ie) 

" D\xTf7\x7\ ' 



623(2:3) = -c 2 [(l + ia)dl - 2 - e»] cos[7r(l - 2te)/4]r(l/2 - «e)^= 



633(2:3) = — ^sign(£ 3 )(l + + a)a^ 

V 71 " Vl x 3| 



(53) 



5.2 The integral representation of the Lazarus-Leblond weight functions. 

The weight function h kp {x, z; x 3 ) (k,p = 1,2,3), corresponding to the work by Lazarus and Leblond 
(1998), are defined as the k-th SIF generated at the point 23 of the crack front by a pair of point forces 
(of opposite direction) exerted on the points [x\ — x, y — , x 3 = z) of the crack faces in the directions 
±e p c| Such a loading can be represented as 

Pi P (x!,x 3 ) = -S ip S(xi - x)S(x 3 - z), (54) 

where Si P is the Kronecker delta. Taking the Fourier transform with respect to 2:3 we obtain 

Pi P (xi,X) = -5 lp 8(xi - x)e lXz , 

and thus, it follows from equation (|51|) . 

hk P (x,z; A) = ie lXz B k j(\)S ip R lh [U J h ] + {-x, A). 

Finally, we deduce 



(55) 
(56) 



h kp (x lZ ;x 3 ) = — I I tB kj (X)6 ip R Ul [ui} + (f3,X)e- l ^e-^-^ x dXd(3. (57) 



00 <j —00 



4 Note that our notation is not identical to the one used by Lazarus-Leblond. In fact hi p (x, 2; x 3 ), h,2 p (x, z; x 3 ) are the 
complex SIF and its conjugate, respectively, h 3p (x, 2; x 3 ) is the mode III SIF. The Lazarus-Leblond's weight functions are 
given by 

, , \ hi p (x, z;x 3 ) + h2 p (x,z;x 3 ) h lp (x, 2; x 3 ) - h 2p (x, 2; x 3 ) 

hj p (x,z;x 3 ) = , hn p (x,z;x 3 ) = — , h U i P (x, z;x 3 ) = h 3p (x, z; x 3 ). 

2 2i 
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The functions hk p {x, z; x 3 ) = hk p (x,x 3 — z) satisfy the parity and homogeneity properties, as outlined 
by Lazarus and Leblond (1998): 

• parity 

— the functions h\ pi hi v (p = 1,2) and h 33 are even with respect to x 3 — z; 

— the functions h\ 3 , /123 and h 3p (p — 1, 2) are odd with respect to 2:3 — z. 

• homogeneity 

— the functions h\ p (p = 1, 2, 3) are positively homogeneous of degree —3/2 — ze; 

— the functions h,2 V {p = 1,2,3) are positively homogeneous of degree —3/2 + ie; 

— the functions h 3p (p = 1,2,3) are positively homogeneous of degree —3/2. 

The constants 7+, 7_, 7m, 7 Z , 7 (see Section are taken from the asymptotics of the weight 
functions hij,hiij,h\iij, as described in Section 3.5 of Lazarus and Leblond (1998). Namely, 



where 
Similarly, 



[hi 2 +ih U2 +i(hn +ih ui )](x,x 3 - z) = \ — \x\ le H + (x, x 3 - z), 



[hj 2 - ih U2 + i(hn - ih ui )](x,x 3 - z) = \ — \x\ le H^(x, x 3 - z), 

Zn 



H + (0,t) = 7+|i| 



{hui2 + ihiui](x,x 3 -z) = \ — \x\ le Hm(x,x 3 - z), 



(58) 



(59) 



[hi 3 +ih U3 \(x,x 3 - z) = \ —H z (x 1 x 3 - z), 



(60) 



where 



him(x, x 3 - z) = \ —H(x, x 3 - z), 
Zn 



H m (0,t) = 7nisign(t)|t|- 2 - 2i£ , H z (0,t) = lz sign(t) |i| ~ 2 " 2 " , H{0,t) = 1 t' 



(61) 



We also note that one can write an equivalent representation of hkp involving the matrix B(x 3 ), 
defined in the text of Section 15.11 



h kp {x,z]x 3 ) = i{B k3 {-) * 5 ip R lh [Ul] + {-x, -)}(x 3 - z) 



(62) 



This gives the exact integral representation of the Lazarus-Leblond's weight functions, which further 
leads to the formulae (fSj) — (JSJ) - The procedure is rather technical and it is outlined in Appendix (|5.3|) for 
one of the constants. 
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5.3 Derivation of the constant 7. 

In this section we will outline the general procedure to obtain the Lazarus-Leblond constants from the 
formula (|57[) . We will derive in detail the formula for the constant 7 defined in (|61|K ; the reasoning 
leading to the other formulae follows. 

We note first that the homogeneity and parity properties of the function H(x 1 x 3 — z) defined in (|61|) -^ 
imply 

H(x,x 3 )= 1 f(-x/\x 3 \)xf 1 (63) 
where /(C) is a smooth monotonic function on R + with the following properties: 

/(0+) = l, /(C) =/oC 2 + 0(C 3 ), asC^+oo, (64) 

where fa is constant. When C < 0, /(C) is defined as zero. 

The property ([64|) ? follows immediately from the fact that the stress intensity factors have to be 
finite as x 3 — > at any x < 0. 

Note that H(x, 0) = jfa/x 2 and 

M*> 0) = J^H(x, 0) = -2=fa\x\- 3 ' 2 . (65) 



Before calculating the constant 7, we will derive the asymptotic behaviour of the function h 33 (x, x 3 ) 



\/\x\/ (2n)H(x, x 3 ) as x — ► . The function h 33 can be written in the equivalent form 

1 \x\ 
h 33 {x,x 3 ) = —=g(x,x 3 ), g(x,x 3 ) = j-=f(-x/\x 3 \)x^ . (66) 
\x\ VZir 



For any fixed x 3 7^ 0, we have 

lim g(x,x 3 )=0. (67) 

On the other hand, for any x < 0, 

g(x,x 3 )dx 3 = -j= [°° f(- x /\x 3 \) l 4dx 3 = 4= f°° f(C)d( = fcf, (68) 

where is a finite constant. Then, (f6"T)l together with (|6"8")l give 

M^a) ~ -^=6(x 3 ), x->Q~, (69) 



which is consistent with the part (d) of the definition at the end of Section [3] 
Using (|57p we can write 



1 r°° poo 1 

h 33 (x ) x 3 ) = - 1[ / -=F3+3(^/|A|)e^ e -^ A dA^, 
4?r z J_ 00 J^oo ^/|A| 



(70) 



where the function 

^ 3 + 3 (£) = -*v1W3i(A)[I7: 3 ] + (£, A) +B 32 (\)[ul 3 ] + (t A) + S 33 (A)[C/! 3 ] + (£, A)} (71) 

does not depend on A, which is easy to check. 

The function -F3 3 (£) is a "+" function and possesses the following asymptotics 

^ + 3 (0 = (i + ^)r 1/2 + ^-^73 3 r 3/2 + o(r 5/2 ), f 6 c+ (72) 
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where 

2 3 (l + e,)- yT^W 
733 = 7= = * ■ (73) 

Now we have to investigate 

/ 2n 

7= lim iJ(x,x 3 )x 3 = lira . / — r/i 33 (x, x 3 )x 3 , (74) 
k->o- x->o- y | a; I 

for any fixed x 3 ^ 0. 

Substituting £ = j3/\X\ in (fTOj) . we obtain 

i />00 />00 

h 33 (x,x 3 ) = — V\Me~ lX3X F+(0e ix ^d^dX. (75) 

J — oo ^ — oo 

Note that the integral in £ is the inverse Fourier transform of the function = /33(£), 

/oo 
fM^dy, (76) 
-oo 

where 

f 33 (y)=0, fory<0, (77) 

and, by the Abelian theorem, 



(78) 



Vv 

Moreover, the function f 33 (y) decays at infinity sufficiently fast. Then (175|) can be written in the form 

i r°° 

h 33 (x,x 3 ) = — ^/\X\e- lX3X f 33 (-x\X\)dX. (79) 

Note that h 33 (x, x 3 ) = for a; > 0, as expected. We also note that the integrand is even with respect to 
A and hence 

1 f 00 

h 33 {x,x 3 ) = - VXcos(\x 3 \X)f 33 {-xX)dX. (80) 



7T JO 

Using the substitution — xX = y > we deduce 



h 33 {x,x 3 ) = 1 , / 7ycos(2/|x 3 /x|)/ 33 (y)dy, x < 0. (81) 



A~x) 3/2 Jo 

Moreover, we can calculate the function /(£) defined in (f63|) as 

H 1 



. ^ cos(y/C) hs{y)dy, C > 0, 
/(C) = 4 V^7C 2 7o (82) 

o, C<o. 

Substituting (JHU) into (JT^J) we obtain 



7 = lim \ — — 

x^o- V 7T x 



Vy/33(2/)cos(y|x 3 /x|)cfy, (83) 
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or, equivalently, for x% separated from zero, we can take the limit as t = |x3/x| — > +00 

j=]imJ-t 2 [ y/yf 33 (y)cos(ty)dy. (84) 
Integrating by parts, we conclude 

7 = lim -\/-* / [%/y/33(y)]'sin(%)dy 

= t Hm -^|[vy/33(2/)]'| J/ =o + ^ 00 [Vy/33(y)]"cos(ty)^| (85) 
= lim n ~ V -IVy/ssCy)]' = 733, 

which completes the required derivation. 

5.4 Comparison of the exact results and asymptotic approximations. 

The exact formulae (J5|)-([9]) are now compared with the asymptotic approximations obtained in the earlier 
work by Lazarus and Leblond (1998b) for small e: 

4u , 8z/log2 2 
8(1 ~ V) i W £-1e + Q(e% 



' tt(2 - v) tt(2 - I/) 

4(1-!/) . 4(l- 1 /)(l-log2) 2 . 
7m = o e + z o e+0(e), (86) 



4 .4(1 + log 2) 2 3 , 
e + i-^r — e 2 + 0(e 3 ), 



7 = 2 ( 2 + ^ +0(e 2 ) 
7 tt(2-z/T K h 



where 



be 



b(b + e) 

It is convenient to use the notation 



(87) 



v = — - — ; (88) 

we note that rj = when the shear moduli of the two elastic media become equal. 

For the sake of illustration, we present the results of comparison for the case when i/+ = v- = 0.3. 
We also consider an extreme situation when v+ = and v— = 0.5. Both diagrams are shown in Fig. [3l 
where we plot the modulus of the ratio of the exact values jj to their asymptotic approximations. 

It is apparent that the asymptotic formulae by Lazarus and Leblond give a very good agreement 
with the exact representations for the case of v + = V- = 0.3. It is also noted that similar accuracy is 
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observed when < v + = v_ < 0.5. The discrepancy becomes more pronounced when v + and f_ are not 
equal, which is illustrated in Fig. [3] for the case of v + = and v_ = 0.5. 

The "worst" constant, in terms of its asymptotic approximation, is 7+. Figure 3] shows the results 
of computations for the same two cases as in Fig. [3] (v + = v- = 0.3 and v+ = 0, V- = 0.5). The 
exact values are shown by the solid line, whereas the dashed line shows the corresponding asymptotic 
approximations . 
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6 Discussion. Application to the wavy crack problem. 



As follows from the earlier work by Willis and Movchan (1995), Lazarus and Leblond (1998) and Bercial- 
Velez et al. (2005), the asymptotic representation for the weight functions can be efficiently used to 
evaluate the perturbation of the stress intensity factors associated with a smooth perturbation of the 
crack front (see the formulae in Section I5TT1 in the text above). In particular, the identity (f2"2"| . written 
in terms of Fourier transforms, takes the form 

a i2 Ri h [ui] + - Ritalin] = -p i R ih [T7i] + , (89) 

where the notations used for components of tractions and displacements are the same as in the text 
above. 

Assume that the crack front is perturbed, within the plane (x\, X3), in such a way that 

x\ < 5<j)(xz), X2 — 0, —00 < 2:3 < 00, (90) 

where S is a small positive parameter, and 4>(x$) is a smooth and bounded function. 

This induces a small perturbation of tractions Aer^ on the plane xi — and associated perturbation 
[Arij] of the displacement jump across the crack, whereas the components of tractions on the crack faces 
can be assumed to remain unchanged. Then the corresponding identity becomes 

{a l2 + A^ a )R lh \U[}+ - RaZZZifik] + [A^,]) = -p^hPl^ ■ (91) 
Subtracting (|89[) from (j9"Tj) we obtain 

A^ i2 Rih - RihZli [Stii] = 0. (92) 

Similar to Section 15.11 one can analyse the delta function term in the left-hand side of (|92| in order to 
obtain the formulae for the stress intensity factors. As outlined in Willis and Movchan (1995), Movchan 
et al. (1998) and Bercial-Velez et al. (2005), the second-order terms in the asymptotic expansions of 
components of the weight functions as well as components of tractions near the crack front are essential 
in this asymptotic analysis. The details of this asymptotic derivation are not the purpose of the present 
paper though - the main result of the present work is the derivation of the closed formulae for the Lazarus- 
Leblond's constants ©-Q. Combined with the asymptotic formulae of Section [2 they complete the 
description of the stress intensity factors near the perturbed edge of the interfacial crack. 
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A APPENDIX. 

A.l Solution of the Wiener-Hopf equation ( 1341) . 

Factorization. We use the factorization obtained by Antipov (1999) for the matrix p~ 1 Go, 

to obtain the factorization of the matrix p~ l G* 

A^G.&signtA)) = F-_i(£,sign(A))r*_(£,sign(A)), 



where 



K + 



1/2 



- fW/2 





At (6 

A.(OAt(0 




£ sign(A) ^ 

'(l + e,)p2(0 (l + e,)p2(0 



sign(A)S+(Q 

sign(A)£+(g) 
5.(£)AJ(0 

_£_ 

^(0 



<5*(£)At(0 

<y.(0At(0 

sign(A) \ 

p2(0 



/ 



Bin(B-(0) sign(A)cos(i?-(0) Ccos(B-(0) 



cos(5-(Q) 
A^(0 



sign(A)sin( J B-(0) £sin(B-(£)) 



P*(0A*(0 



P*(£)A*(£) 



(A.1) 



(A.2) 



(A.3) 



(A.4) 
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where 



^(0 



(A.5) 



7T/9*(£) 



1 



47T 



log 



^(t 2 + l) 1 /2 + l 



dt 



d»(t 2 + l) 1 /2_l( < 2 + 1 )l/2( t _^- 

and the contour of integration is as follow: 

r = U = a - *0, -oo < a < -|o|} U U = a + «0, |a| < & < oo} U {£ e R, |£| < |a|}. 
Other functions involved in the matrices (| A.3|) — (|A.4p are defined by 

E+(0 =id»[cosB+-id»p»sinB+], E+ 2 (f) = id.[— sinS+ + id* cosB+], 



(A.6) 



where the constants used here and in the sequel are 



d = y/l - d 2 , di = log - 



1 + d* 



d* ' d V 2 ■ 



(A.7) 



The branch cut for the logarithmic function in (|A.5[) a is defined in the same way as in Antipov (1999), 
formula (3.9). 

The inverse of matrices (|A.3[) and (IA.4|) are as follow 



/ 



Y-l = (^ + i)-^ 



-(l + e.)£ 



(l + e*)sign(A) ^A+E+ a 



dl 

Si § n ( A ) 



sign(A) 



A+V+ 



dip: 



(A.8) 



/ 



^ = (e-*) 1/2 



-/9*(0Ar(0sinB7(0 A;(OcosB-(0 \ 



£ sign(A)A-(OcosiJ-(0 ^#A;(£) sini?-(0 



sign(A) ^A-(OcosB-(0 



V 



A-(OsinB-(0 



(A.9) 



/ 



r/te solution. Substituting the representation for the matrix G*(£, sign(A)) into the eq. (I34|) . we 
obtain 

^.+(00^(0 = ^.-(0^(0- (A.io) 



21 



Let us introduce the notations 

e* + = y»+(O0+(O, £*- = y-(0^(0- 



(A.ll) 



It follows from the asymptotics of the matrices V*±(£) at the points ±a and ±i and from the asymptotics 
of the weight functions at infinity that 

/ ^ + ^L \ 



C3 C4 



C 5 



(A.12) 



Analyticity of the vectors (£) at the points ±i and ±a requires the following conditions to be 
satisfied 



(1 + e»)sign(A)Ci 



C 5 



= 0, 



(A.13) 
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C* 2 + i sign(A)M_C 3 = 0, iC 4 - d*C 5 = 0. 

Note that the three linearly independent weight functions ([5r?)) - P0)) have been obtained from (|A.12[) , 
using (|A.13[) and the following choice of the coefficients Cj : 

(1) for the first weight function: 

Ci = (1 + e*)-\ C 2 = C 3 = C 4 = C 5 = 0, C 6 = -Ud*M- sign(A), 

(2) for the second weight function: 

d = 0, C 2 = -id*ik£ sign(A), C 3 = d.Af_, C 4 = C 5 = C 6 = 0, 

(3) for the third weight function: 

Ci = C* 2 = C 3 = 0, C 4 = -»M_, C 5 = M_d-\ C 6 = -d- x M_(i - a) -1 . 

A. 2 The asymptotics of weight functions. 

To obtain the asymptotics of the weight functions, as £ — > 00, we will make use of the following repre- 
sentations 



2£ , ^ /log? 



7r£ ±Z 



^(0 = I 1 log? + log A + Bo + O 
Zir Zir ±1 



cos Bf (0=^(0+0 



log? 



log? 



(A.14) 



sinBf(C) =^(0+0 



logC 



where 



So 



1 

2Wo 
1 



log 



d*(t 2 + l)V 2 + i 



(t 2 + l)Va 2 



i , 1 — ad* 
log — j 



(A.15) 



^(o = §(4^r + ef^r*), ^(0 = -2(4 - efa-^-*). 
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We can show that 



so that 



1 r c 

2Wo 



log 



d,(t 2 + l)V2 + l 



d*(t 2 + l) 1 /2-l 



di 



(*2 + 1)1/2 4' 



In this case 



„ 7T i 1 — ad* 



(A.16) 



(A.17) 



where 



*io(0 

^20 (0 

^20 (0 



'1 


— za 


2 


/l 


— ia 



(-«U + (0)> F n (0- 
(id*c + (£)), F 21 (0 = 



1 — ia 



c + (0--d.(i-a)8 + (Z) 



1 — ia 



s + (0+2 rf *(*- ft ) c+ (0 



r/ 2 
"() 



d 2 



/l 


— ia 


2 


'1 


— za 



„ 1 / 1 — ia I i — a 

*-(0, P n(0-^ 



(A.18) 



2 C «>' P2l(?) = d 2 



and 



1 — ad* 



(A.19) 



Finally, we can write the leading terms in asymptotics of the Fourier transform of the weight functions, 
which are of the order 0(^ -1 / 2 ) and 0(^ 1 / 2 ), as £ — > 00, for the displacement and traction fields, 
respectively. Namely 

(1) The first weight function: 

p\ x ] ~ sign(A)r 1/2 {-2fio(0 + l [1 - 2Fn(0 + i*io(0]} , 



[C/l 2 ] ~ sign(A)r 1/2 {2^20(0 + j [2F 2 i(0 - ^20(0]} , 



(A.20) 



-1/2 



2P 10 (£) 



J *12 S1 g 



n(A)£ 1/2 (2id» 



1 + e* 



+ 2id*Pn(0+d*Pi (0 



-sign(A)e 1/2 |2id*P 20 (O + 1 [2id*P 21 (0 + d*P 20 (0]} , 



(A.21) 



J *32 



^1/2 



1 1 

1 + e* + I 



2(1 + e,) 



-2id*Pi (£) 
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(2) The second weight function: 

Pli\ ~ r 1/2 {^20(0 + 1 
[^ 2 ]~r 1/2 (^io(o + 7 



+ iJ21 (0 + ^20(0 



2d 2 l d*(i - a) 



^11(0+2^(0 



[C7 2 3 ] ~ sign(A)r 1/2 



1 + e* 



2dod*(i-a) £ 



, M 3. 1 + e* 
**2o(0 + 



2 2d 2 ,d*(i - a) 



~ o /2 \d*p 20 (o + \ 



^7rA-. r + d*P 21 (0 - lid*P 20 (0 

2dgd*(l — a) 2 



v 2 

^*22 



E 2 32 ~ sign(A)0 /2 



Ao(0 + i 



-rf*Pii(0+ 2 jrf * p i"(0 



} 



1 

+ T 



2d 2 d* (i - a) £ 



3i 1 1 

4dgd*(i - a) J J 



(3) The third weight function: 



*2J 



, r 1/2 F2o(o + — — ,f w (o 

d* [ d*(i-a) 



*2l(0 + 



-(«-J)j*o(fl + (§+**, i( ._ a) 



d*(i - a) 
1 



^ii (0 

fio(0 



j-r 1/2 f fi (o - , * . f 20 (o + 



-("-2)^10(0 



+n (0 - 

3 



d*(i — a) 
1 



„ + la , 

2 / a* (i — a) 



F21 (0 

^20 (0 



sign(A) ^_ 1/2 f 1 
d* \£ 



^20 (0 



d*(« — a) 



fio (0 



J *12 



^*22 



o /2 -ifl»(o 



Ao(0 , 1 



d*(i — a) £ 



e i/ 2 | 



^io(0 



S: 32 ~ sign(A)0 /2 ^ - 



^21(0-^+^^0(0 

; ^11(0-^^^0(0 



d*(i — a) 



2d*(« - a) 



^20 (0 , 1 



d* (i — a) £ 



-^20(0 



1 



^11(0+^2+^^0(0 

:i ^(0-^-^^20(0 



d*(i — a) 



2d*(z — a) 



d* (« — a) 



Ao(0 
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A. 3 Asymptotics of stress field. 

Stress components have the following asymptotics: 



i f T r, n -1/2+ie \ -1/2-ie . ,, n 1/2+ie , * / s 1/2— 1 

cri2 — 7 =\K(x 3 )x 1 - K (x 3 )x 1 + A(x 3 )x 1 l - A (x 3 )x 1 / \, 

ZV Z7T I- J 

a 22 ~ — ^= {^(x 3 )x- 1/2+l£ + iH^K 1 / 2 ^ + Afo)*}' 2 ^ + A*^)^ 2 ^} , (A.26) 
2\/27r l- J 



2V2tt 

Km(x 3 ) -1/2 , Am(x 3 ) 1/2 
cr 32 7=^1 H ^^ x l 1 



/2tt V2tt 
The corresponding Fourier transforms are: 



ai 2 ~ i (^(A)-i-r 1/2 "" - iT (A)^/T 1/2+J£ + A(A)-^-/3- 3/2 - ie - A*(A)^/r 3 / 2 + 
2 t 2cie 2c 2 251 e 2# 2 

l - (^(A)-^/3- 1/2 - 4e + K*Wp-0~ 1/2+l£ + A(X)-^—r 3/2 - ie + A*(X)^r 3/2+l£ \ , 
I y 2cie lc 2 zgi e 2g 2 J 
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*32 ~ ^Kui(\)r 1/2 + ^^A m (\)f3- 3/2 , 
where 

1 2r(l/2 + ie)' 2 2r(l/2-ie)' ° 
(l-i)v^ (1 - *)V5F 



(A.27) 



(A.28) 



2r(3/2 + ie)' a " 2r(3/2-ie) 

A. 4 Asymptotics of displacement components. 

The crack opening displacement are characterised by the following asymptotics 

L J ^Fcosh(Tre) \l + 2ze V U 1 - 2ie 1 17 

+ B^X-^) 3 /^ - B*( a : 3 )(- a;i ) 3 / 2 -"| , 



+ fifoX-si) 3 / 24 * + B*(x 3 )(- a;i ) 3 / 2 -"| 



[«s] - {^iiiM(-xi) 1/2 + B m (x 3 )(- Xl )^} . 
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The corresponding Fourier transforms are: 

[ui] ~ -i lK(\) Vl e^/ 2 -"- - K*(\)^f3- 3/2+te + B(X)w ieo (3- 5 / 2 -^ - B* (X)—/3~ 5 / 2+te X 
I eo e J 

[u 2 ] ~ K (A)wie 0/ 9- 3/2 ~ 4£ + K*(\)—f3- 3/2+l€ + B{X) Wl e l3- 5 / 2 - l£ + B*(X)— /r 5/2+ ", 

e e 



[«s] ~ 2(6 + e) j-i±^ m (A)/r 3 / 2 + ^2±is m (A)/3- 5/2 } , 



(A.30) 
where 

ibS ibdZ 

4ci 4C2 (A.31) 



wi = ^ r(3 + 2ie)j W2 = ^ r(3 _ 2z£) , 



2G 



